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Switching of Nonequilibrium Depletion Force Caused by Blockade Effect 
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The concentration-dependent switching of the non-equilibrium depletion forces between obstacles 
in an interacting Brownian gas flow is presented. It is shown that this switching is caused by 
the blockade effect for the gas particles. With increasing equilibrium gas concentration, the gas 
particles blockade causes the obstacle wake inversion (trace profile “turn-over”) that, in turn, leads 
to the change of sign of dissipative interaction. Some non-linear effects such as formation of a 
cavity-like sparse wake behind the obstacle and the dissipative pairing effect are discussed briefly. 
The results are obtained within the lattice gas model in the mean-field approximation. 


I. INTRODUCTION 

The motion of impurities in a Brownian gas is accom¬ 
panied by its density perturbation (impurity trace, or 
wake) that, in turn, induces additional non-equilibrium 
correlations between the impurities. These correlations, 
called dissipative or non-equilibrium depletion forces [1 
4], are responsible for the coherent part of the collective 
friction force as well as for possible dissipative structures 
of impurity ensemble. 

In equilibrium, the depletion interaction is usually 
short-range; the spatial range is of the order of the deple¬ 
tion agent characteristic length scale [5, 6]. In contrast, 
the non-equilibrium forces between the impurities may 
exhibit long-range behavior due to a long-living diffusion 
trace induced by their motion [7-11]. In addition, such 
forces often have unusual properties, e.g., they violate 
the Newton’s third law [1, 2, 12, 13]. 

It is clear that specific properties of the dissipative in¬ 
teraction are determined by the structure and the shape 
of the perturbation profile, which, in turn, significantly 
depends on the interaction between the gas particles. 

The presence of inter-particle interaction, even short- 
range one (as in a lattice gas, when a lattice site can be 
occupied by only one particle), is responsible for a num¬ 
ber of unexpected kinetic effects, e.g., “back correlations” 
effect [14], drifting spatial structures [15-17], effects of 
“negative” mass transport [18-20], dissipative pairing ef¬ 
fect for tracers passing through a lattice gas [21]. 

An important mechanism often underlying such kinetic 
phenomena is a blockade effect. As was shown in [22], 
the blockade effect implicates significant changes in the 
form of an obstacle trace. The trace structures for the 
cases of impurity moving in a lattice gas and fixed impu¬ 
rity or obstacle in a gas flow do not coincide in general. 
In particular, the obstacle trace can take atypical form 
characterized by switching of its direction depending on 
the equilibrium gas concentration. Usually, the profile of 
the gas perturbation induced by moving inclusion has a 
dense localized region ahead of the inclusion and a long 
low-density tail behind it, which damps according to a 
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power law, e.g., |7 11]. However, for an obstacle in a gas 
flow, the gas perturbation profile can take reversed struc¬ 
ture with an extended dense region ahead of the obstacle 
and a localized low-density region in its wake [22], It 
should be noted that similar effect in one-dimension was 
discussed in [23]. This switching of the profile directivity 
is caused by the effect of particles blockade ahead of the 
obstacle. The blockade is due to hard-core-like interac¬ 
tion between the particles and is possible if the flowing 
gas density is high enough. 

As the inversion of density perturbation takes place 
for arbitrary number and configuration of obstacles, see 
[22], we may expect that this inversion has to affect the 
dissipative interaction between the obstacles. 

The aim of this paper is to show that trace inver¬ 
sion is responsible for the change in sign of the induced 
non-equilibrium interaction between the obstacles, e.g., 
switching of the effective attraction to repulsion. 

We consider the simplest model of a lattice gas with 
two obstacles whose sizes are much larger than the lattice 
constant. We use the mean-field and the long-wavelength 
approximations, neglecting short-range correlations and 
fluctuations in the gas, see [22], 


II. TRACE INVERSION IN THE LATTICE GAS 
MODEL 

As was shown in [18, 22], the problem of an obstacle 
in a lattice gas flow can be considered as a limiting case 
of a two-component gas, one of the components is at rest 
while another one is flowing in a uniform external field. 

We resort to the simplest model of a two-component 
lattice gas, when each lattice site can be occupied only by 
one particle, see [14]. The kinetics of multi-component 
lattice gas is defined by the jumps of particles to the 
neighboring vacant sites. The variation of the i-th site 
occupancy by the particles of a sort during the time in¬ 
terval At, t 0 <C At -C t i (r 0 is the duration of a particle 
jump to a neighboring site, r; being the lifetime of a par¬ 
ticle on a site), is described by the standard continuity 
equation (see, e.g., [14, 24]) 

<(t + At) - <(t) = £ (J“ - J“) + 6J?, (1) 
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where nf = 0,1 are the local occupation numbers of par¬ 
ticles a at the i- th site, ^1 — Yp n j ) At gives 

the mean number of jumps (of a particles from site i to 
a neighboring site j per time At, /3(a) indexes particle 
species), = v a is the mean frequency of these jumps. 
The term 3.7“ = Yj($Jji ~ £J,“) stands for the Langevin 
source that is defined by the fluctuations <5 Jj* of the num¬ 
ber of jumps between sites j and i during At [24]. These 
fluctuations are caused by the fast, as compared to the 
time scale At, processes and will be neglected for simplic¬ 
ity. It means that we neglect fluctuation-induced forces 
that may be appreciable for closely located inclusions. 

In what follows we consider only two components, de¬ 
noted by n and v. In the absence of external fields for 
the regular lattice we suppose that idt = v = const for 
the mobile component n, while the component v is as¬ 
sumed to be at rest, z// = 0. The presence of driv¬ 
ing field leads to asymmetry of particle jumps. Suppos¬ 
ing the activation mechanism of jumps and the driving 
field G to be weak, the frequency may be written as 
z/tt « z/[l + (G, r i~ rj)/(2kT)\, or i/ ± m i>±8v, where v + 
denotes the jump frequency along the field, v~ - against 
it, 5v = va\Q\/(2kT) (a is the lattice constant), condi¬ 
tion a|G|/(2fcT) < 1 being satisfied. 

Equations for the average local occupation numbers 
can be obtained from Eqs. (1) using the local equilibrium 
approximation (Zubarev approach) [24, 25] which coin¬ 
cides, in our case, with the mean field approximation [26]. 
Introducing time derivatives [27], in the long-wavelength 
approximation, see [15-18], macroscopic kinetics of the 
mobile component n is given by the equation 

d T n = V 2 n — V(uVn — nVu) — (g, V)[n(l — v — n)], (2) 

where n = n(r,r) and v = v(r) are the average occu¬ 
pation numbers of the two components at the point r, 
g = aG/(2kT). Here, we have introduced dimensionless 
spatial coordinate r/o —> r and time r = vt, and d T 
stands for the partial time derivative. 

Second term on the right-hand side of Eq. (2) is a re¬ 
duced nonlinear mixing flow which describes mutual drag 
of one component by another and arises as a result of the 
particle distinguishability and the local repulsion (due to 
excluded volume constraint) in a lattice gas [18|. This 
flow leads to the array of anomalous diffusive transport 
effects, e.g., the drag effect, formation of the drifting spa¬ 
tial structures [15-17], effects of “negative” mass trans¬ 
port [18-20], and induced long-time correlations [28]. 

In order to demonstrate the trace inversion effect men¬ 
tioned above we now consider a single macroscopic ob¬ 
stacle of a circular form in a flow of driven lattice gas 
particles. Fig. 1 presents the two-dimensional numerical 
solutions of Eq. (1) in the mean-field approximation for 
different equilibrium concentrations no- 

For low concentrations (no < 1/2), the structure of 
the profile is typical [1, 3, 4, 11] and is characterized by a 
lengthy depleted region behind the inclusion (wake) and 
a localized dense region in front of it, Fig. la. 
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Figure 1: The two-dimensional equilibrated concentration dis¬ 
tributions Sn(x, y) = n(x, y) — no (mean occupation numbers) 
of the lattice gas particles near the inclusion at the different 
values of equilibrium concentration no are presented on con¬ 
tour plots: (a) no = 0.2, (b) no = 0.5, (c) no = 0.8. The exter¬ 
nal field g (|g| = 0.5) is directed along the i-axis; the imper¬ 
meable (no = 1) circular inclusions of radius R = 7 (in units 
of a) are placed at the origin (denoted by circles). The gray 
background corresponds to the equilibrium gas concentration 
no for every contour plot, in consistence with the colorbar. 
Dashed contour lines underline the exact inversion between 
profiles (a) and (c), while the corresponding one-dimensional 
profiles 6n(x) = n(x, y = 0) — no show the “switching” of 
asymptotic behavior. 



At high concentrations (no > 1/2), the profile takes 
unconventional form with an extended dense region in 
front of the inclusion and a localized depleted one in its 
wake, resembling the form of a cavity, Fig. lc or Fig. 2. 
Hence, the main part of the density perturbation is now 
shifted to the region in front of the obstacle and directed 
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upstream. 

This behavior is caused by the blockade effect of the gas 
particles flow due to the short-range repulsion between 
them (the condition that a lattice site can be occupied 
by only one particle). For sufficiently high concentrations 
no, gas particles have no time to leave the blockade zone 
ahead of the obstacle via lateral diffusion and, as a result, 
the dense region ahead of the obstacle has to grow. The 
blockade effect is non-linear and becomes significant near 
the obstacle surface and/or for large obstacles whose size 
is much larger than the lattice constant. 

The trace inversion directly follows from the symmetry 
of Eq. (2) that can be represented in the form 

d T n = V(/iVn — nVh — gn/i), (3) 

where h = 1 — m — v is the vacancy concentration, i.e., 
concentration of the empty lattice sites. Equation (3) is 
invariant under transformation n h, g —g and can 
easily be rewritten as 

d T h = V(nVh — hVn + gnh). (4) 

Equation (3) describes the kinetics of the gas particles 
in sweeping field g at equilibrium gas concentration n o 
while Eq. (4) corresponds to the vacancy kinetics in the 
opposite field —g at the equilibrium concentration of va¬ 
cancies ho = 1 — no- The tail of the obstacle trace in the 
gas flow is characterized by a depleted gas region behind 
the obstacle, whereas the trace in the corresponding va¬ 
cancy flow (which has opposite direction) has a depleted 
tail of vacancies that corresponds to a dense region of gas 
particles. 

For distribution v(r) with mirror symmetry along the 
field direction r|| || g, i.e., ■u(r||) = v(— r||), each density 
perturbation profile 8n{y,no) = n(r, no) — no at a given 
concentration no can be obtained by the inversion trans¬ 
formation of the profile 8n( r, 1 — no) at concentration 
1 — no, that is expressed by relation 


M r || + rj_; 1 - n 0 ) = -Sn{- ry + rj_; n 0 ), (5) 

r || and r± are longitudinal and transverse components of 
r, with respect to the external field g. 

For relatively small obstacle with size R, the asymp¬ 
totic behavior of the gas density perturbation 6n(r) at 
|r| R as well as the trace inversion effect can be 
obtained in the linear approximation, see [22]. Far 
from the obstacle, one may assume that the distribution 
n = no+Sn should weakly differ from the equilibrium one 
n 0 . Thus, it is possible to linearize Eq. (3) by neglecting 
term with ( Sn ) 2 . For the case of two-dimensional lattice 
gas the trace asymptotic Sn(r) for a single obstacle of 
radius R is given by (see [22]) 


x ^ ( (X n 
n=0 
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Figure 2: The example stationary wake (numerical result) 
whose profile m(x, y = 0) is of the form of Bloch wall, describ¬ 
ing the cavity behind the inclusion (black circle); equilibrium 
medium concentration mo = 0.63, field |g| = 0.5, inclusion 
radius R — 7. 


where $ is the angle between g and r, r = |r|; constants 
a n = a n {R) are determined from boundary conditions 
on the obstacle surface, /3 = (1/2 — n 0 )/|l/2 — no| = ±1, 
and q = |l/2 - n 0 ||g|. 

In the particular case of no = 1/2, Sn is an odd func¬ 
tion of x, see (5), and the density perturbation becomes 
asymmetrical and extended in the ±x directions, see 
Fig. lb. As was shown in [22], its asymptotic behavior 
damps according to the power law 6n( r) ~ —(g, r )R/r 2 
that is similar to the scattered electrostatic potential for 
a dielectric particle in a uniform electric field. 

The linear flow approximation gives the qualitative de¬ 
scription of the trace inversion effect and the character 
of the asymptotic behavior of Sn. In the particular case 
of the point-like inclusion the method gives Sn ~ r -3 / 2 
that is in satisfactory agreement with numerical results 
and coincides with asymptotic behavior of the relaxation 
of the trace induced by a moving intruder presented in 
[9], 

However, as was mentioned above, non-linear effects 
may be significant near the obstacle surface. Fig. 2 shows 
the role of this non-linear effect in the formation of the 
density perturbation profile Sn. The behavior Sn near 
the surface has a pronounced step-like character that de¬ 
scribes a localized dense region ahead of the obstacle at 
no < 1/2 or a localized depleted region behind it (resem¬ 
bling the form of a cavity) at no > 1/2. 

Intuitively, such behavior can be expected. In the par¬ 
ticular case of a lattice gas with no obstacle Eq. (3), with 
v = 0, admits the kink-like solution. 
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III. NON-EQUILIBRIUM FORCE IN A 
LATTICE GAS MODEL 

In the non-equilibrium case, there are various ways to 
introduce dissipative force or interaction between inclu¬ 
sions via Brownian gas environment. The ways are not 
equivalent to each other and may lead to different results 
in general, see [2], In particular, there is the Smolu- 
chowski approach, the approaches based on a master 
equation, free energy, etc., see [1, 2, 29, 30]. In this pa¬ 
per, we use the approach based on a master equation, 
see [2], that is close to the often exploited method of fast 
variable elimination, see, e.g., [31 33]. Considering inclu¬ 
sions as particles of a heavy gas component with given 
interaction potential between them and the original gas 
particles (the light component), the force (effective inter¬ 
action potential) between the inclusions can be obtained 
by averaging it over the light component. 

In the equilibrium case, the average force acting on 
the jth inclusion centered at the lattice site Rj can be 
formally written as (see [2] for details) 



= -E< n *>V% y (r*-R 3 ) (7) 

i 

where rii = 0,1 is the occupation number of site r,. The 
jth inclusion is modeled by the interaction potential V 
between its center and the gas particles. 

(rii) = ^ {nip({n},°)} (8) 

{nil 

is the averaged occupation number (gas concentration) of 
the ith site and p({n},0) is the equilibrium probability, 
or statistical operator in matrix representation [24], of 
finding a given occupancy configuration {n} 

p({n}, 0) = Z~ x exp (~H{n}/kT), (9) 

Z = ]T]{ n } exp(— H{n}/kT), where H{n} is the lattice 
gas-inclusions Hamiltonian 

H{n} = H 0 {n}+ ^2V(ri (10) 

id 

and Hq{u} is the lattice gas Hamiltonian without inclu¬ 
sions. 

In the continuum limit and for inclusions given by a 
hard-core-type potential, the expression for the average 
force takes the custom form (see, e.g., [29]) 

F f = - J s n(r)(n(r)) dr, (11) 

where Sj is the jth inclusion surface and n(r) is its exte¬ 
rior normal at the point r. Expression (7) gives standard 


form for the equilibrium depletion force that can be ex¬ 
pressed in terms of the free energy of a lattice gas with 
the volume occupied by impurities being excluded [1, 2]. 

In the non-equilibrium case the dissipative or non¬ 
equilibrium force can be written in a form similar to that 
of (7), see [1, 2] 

*j(t) = -Y l [Vn. s V(T i -R j )] 5{m) t 

i 

-A — / n(r)S(n(r)) t dr, (12) 

Js 

where S(rii)t = (rii)t — (rii) is the density perturbation of 
the lattice gas and 

(ni)t = '^2{n z p{{n},t)}. (13) 

{«} 

p({n},t) is the probability of finding a given occupancy 
configuration {n} at the time t, which satisfies the master 
equation for hopping process (see, e.g., [34]) with initial 
conditions corresponding to the equilibrium distribution 
(9). 

The total dissipative force, Eq. (12), exerted on an 
inclusion is defined by the induced density perturbation 
(obstacle trace) and can be conventionally represented by 
two parts. The first part is individual friction force on 
the isolated inclusion. The second part takes into account 
collective, or coherent, contribution that is related to the 
dissipative interaction between the obstacles, caused by 
their mutual influence via their traces. 

In what follows, we will be interested in the case of 
steady-state non-equilibrium interaction, i.e., the limit¬ 
ing case t —> oo. 

To describe the time behavior of averaged occupation 
number ( m)t we exploit an approach based on the kinetic 
equation for dynamical variable n,, Eq.(l), and use the 
results for the average occupation number n obtained in 
the previous section. 


IV. FORCE SWITCHING 

In this section we will show the change of sign of dissi¬ 
pative interaction between two obstacles, that is caused 
by their trace inversion. First, we qualitatively show that 
such interaction switching with increasing gas concentra¬ 
tion n o is a direct consequence of the system symmetry 
(5). To demonstrate this property we resort to the case of 
two inclusions located perpendicular to the gas flow, see 
Fig. 3. For simplicity, the inclusions are given by smooth 
distributions with compact carrier Vj( r) = v(r — Rj), 
the inclusions’ centers are located at points Ri = ( 0,yo ) 
and R 2 = (0,—yo). In this case, the force exerted 
on the jth inclusion by gas perturbation Sn( r) reads 
F j = J 5nX7vj dr. To show the possibility of the force sign 
changing, it is enough to consider only the y-component 
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Figure 3: Transverse alignment. The particle concentration distributions at three different regimes of the dissipative interaction: 
(a) no = 0.2 — effective repulsion (F-/ 2 > 0, F 21 < 0, \F 21 \ = |F/ 2 |), (b) no = 0.5 — no interaction (F 21 = F^ 2 = 0), (c) no = 0.8 
— effective attraction (F 21 > F/ 2 ). The external field g (|g| = 0.5) is directed along the a;-axis; the impermeable (no = 1) 
circular inclusions are of the radius R = 7 (in units of a); the distance between inclusion centers d = 2yo = 4 R. The gray 
background corresponds to the equilibrium gas concentration no for every contour plot, in consistence with the colorbars. 


of the force Fj 7 (no) exerted on an inclusion at equilib¬ 
rium gas concentration no- Taking into account (5) one 
can obtain 

Fj(n 0 ) = J 5n(r;n 0 )d y Vj(r)dr = -F/(l - n 0 ). (14) 

In the particular case of half-filled medium (no = 1/2) 
F?( 1/2) = —Fj( 1/2) = 0, i.e., the effective interaction 
between the obstacles vanishes. 

We next consider numerically the wake-mediated force 
between two obstacles for two configurations: the ob¬ 
stacles are placed parallel and perpendicular to the gas 
flow. The total force exerted on a given obstacle includes 
the part associated with individual friction force and that 
associated with the influence of another obstacle. To sep¬ 
arate out the inter-obstacle contribution from the total 
dissipative force we consider the quantity fl] 

Fy = F; - F° = f [tfn(r, R,, R,) - Sn( r, RQ] Vu i (r)dr, 

(15) 

where F* is the total force acting on the ith obstacle in 
the presence of jth one and F/ is its individual friction 
force. 

Transverse alignment (Fig. 3). One can see from 
the symmetry of this configuration that the forces two 
obstacles exert on each other are equal and opposite, 
F/ 2 = —F|i, i.e., their interaction is always Newtonian, 
see Fig. 4. 

At low equilibrium concentrations (n 0 < 1/2), Fig. 3a, 
the dissipative interaction manifests itself as an effec¬ 
tive repulsion between the obstacles, since F 2l < 0 and 
Ff 2 > 0, see Fig. 4. Qualitatively, this repulsion is simply 
explained by the overlap of density coats around the ob¬ 
stacles that leads to the onset of a dense region between 
them acting like a repelling barrier , see Fig. 3a. 

In contrast, at no > 1/2, the overlap of the individual 
density perturbation coats of the obstacles results in the 
formation of the extended dense zone ahead of them that 
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Figure 4: The dissipative forces F/ 2 (rio) and F 21 (no) as a 
functions of the equilibrium medium concentration no in the 
case of transverse alignment. The values marked with circles 
correspond to the profiles presented on Fig. 3. 


blocks the gas flow, so that the region between the obsta¬ 
cles becomes depleted. As Fig. 4 suggests, this collective 
blockade effect of gas particles leads to the effective at¬ 
traction between obstacles in a dense medium, F 2 \ > 0 
and Ff 2 < 0. 

Thus, when the gas concentration no increases the dis¬ 
sipative interaction between the obstacles switches from 
repulsion to attraction. In addition, the non-linear inter¬ 
obstacle attraction may characterize the pairing effect ac¬ 
companied by the creation of common perturbation coat 
around the obstacles. The effect of a similar nature was 
obtained earlier in [ 11] for two driven intruders. 

In the no = 1/2 case the effective interaction between 
the inclusions vanishes, Ff 2 = F 2l = 0, irrespective of 
the inter-inclusion distance. 

Longitudinal alignment (Fig. 5). In the case of low con¬ 
centration (no < 1/2), the typical situation for Brownian 
systems is shown. An inclusion falling on the depleted 
wake induced by another inclusion is effectively attracted 
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Figure 5: Longitudinal configuration. The lattice gas parti¬ 
cle concentration distributions n(x, y) corresponding to three 
different regimes of the dissipative interaction: (a) no = 0.2 
— effective attraction (|Ffj| > |Ff 2 |), (b) no = 0.5 — anti- 
Newtonian interaction (F|j = Fj 2 ), (c) no = 0.8 — effective 
repulsion (|Ffi| < |Ff 2 |). The external field g (|g| = 0.5) is 
directed along the x-axis; the impermeable (no = 1) circular 
inclusions are of the radius R = 7 (in units of a), their posi¬ 
tions being marked with black circles; the distance between 
inclusion centers d = 2xo = 10F. The gray background cor¬ 
responds to the equilibrium gas concentration no for every 
contour plot, in consistence with the colorbar. 


to it since the friction force in the more depleted region is 
weaker [1, 3]. This type of effective interaction is referred 
to as the wake-mediated [35]. As Fig. 6b suggests, the 
second obstacle does not practically affects the first one, 
Ff 2 ~ 0. In contrast, at high concentrations no > 1/2, 
the second obstacle does not feel the influence of the first 
one, ss 0, whereas the first obstacle comes under 
the excess pressure of the dense gas region created ahead 
of the second one due to the blockade effect. As a re¬ 
sult, the effective interaction changes its sign switching 
from attraction to repulsion, Fig. 6b. In particular case 
of n o = 1/2, the effective interaction between the inclu¬ 
sions becomes strictly anti-Newtonian, Ff 2 = F/j ^ 0 
(see Fig. 6b), which holds for arbitrary inter-inclusion 
distance. 

Note that for a dense gas in the blockade regime, the 
second obstacle “pushes” the first one upstream, thus re¬ 
ducing the total friction force Ff exerted on the first 
obstacle, Fig. 6a. 

The dissipative interaction between the inclusions nat¬ 
urally vanishes in the limit of strongly sparse (empty) 
medium no —t 0, due to wake depletion. The same is 
true in the total jamming limit n o —t 1. 


Figure 6: (a) Concentration dependence of the dissipative 
forces Fq (no), Ff(no), and F 2 (no), (b) The forces Ff 2 (no) 
and F|j (no), acting between the inclusions in the longitudinal 
configuration. The values marked with circles correspond to 
the profiles presented on Fig. 5. 


V. CONCLUSIONS 

In this paper we have shown the switching of the non¬ 
equilibrium depletion force between fixed impurities or 
obstacles immersed in a flowing and interacting Brownian 
gas. To this end we have considered the simplest model 
of a lattice gas with two obstacles whose sizes are much 
larger than the lattice constant. We also used both the 
mean-field and the long-wavelength approximations, ne¬ 
glecting short-range correlations and fluctuations in the 
gas. 

The non-equilibrium interaction between the obstacles 
is caused by the density gas perturbation or wakes in¬ 
duced around them by the gas flow. The force switching 
is due to the blockade effect of the gas particles and man¬ 
ifests itself by changing the direction of force to its oppo¬ 
site. With increasing the equilibrium gas concentration, 
the blockade effect provokes the wake inversion (or the 
wake profile turn-over) that, in turn, leads to switching 
of non-equilibrium depletion interaction, e.g., from the 
effective repulsion to attraction. 

In contrast to the equilibrium case, the non¬ 
equilibrium depletion force exhibits long-range character 
due to extended wake tail behind the obstacles which 
damps according to a power law which can be estimated 
in the linear approximation. However, non-linear effects 
become significant near obstacle. The density pertur¬ 
bation near its surface exhibits a pronounced step-like 
profile that describes the formation of a sparse cavity-like 
region of the gas behind or a dense drop-like region ahead 
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the obstacle in the cases of high and low equilibrium gas 
concentration, correspondingly. In turn, this can lead to 
the effect of dissipative pairing for two closely located 
obstacles that accompanied by the creation of a common 
coat of gas density perturbation around them. However, 
with using the mean-field approximation we lose infor¬ 
mation on short-range correlations which may be signifi¬ 
cant near the obstacle surfaces and in the case of closely 
located obstacles, see [14, 21]. Moreover, neglecting fluc¬ 
tuations in a gas, e.g., the term in Eq. (1), we do not 


take into account the fluctuation-induced (Casimir-like) 
forces, see, e.g., [29, 36—41). 
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